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We study properties of two different atoms at a single optical lattice site at a heteronuclear atomic Feshbach 
resonance. We calculate the energy spectrum, the efficiency of rf association and the lifetime as a function of 
magnetic field and compare the results with the experimental data obtained for 40 K and 87 Rb [C. Ospelkaus et 
al, Phys. Rev. Lett. 97, 120402 (2006)]. We treat the interaction in terms of a regularized S function pseu- 
dopotential and consider the general case of particles with different trap frequencies, where the usual approach 
of separating center-of-mass and relative motion fails. We develop an exact diagonalization approach to the 
coupling between center-of-mass and relative motion and numerically determine the spectrum of the system. 
At the same time, our approach allows us to treat the anharmonicity of the lattice potential exactly. Within the 
pseudopotential model, the center of the Feshbach resonance can be precisely determined from the experimental 
data. 

PACS numbers: 34.20.Cf, 34.50.-s, 37.10.De, 03.75.Kk 



INTRODUCTION 

Motivated by the intriguing perspectives of heteronuclear 
molecule formation, observation of charge-density waves |Q]], 
boson-induced fermionic superfluidity J2, d, 01 in optical lat- 
tices (H, and supersolids [6], Fermi-Bose mixtures have re- 
cently attracted lots of attention. An important step in this di- 
rection was the simultaneous trapping of bosons and fermions 
in a three-dimensional (3D) optical lattice 10, [Hi- Recently, 
even heteronuclear molecules IH (T^l were created by means 
of a magnetic field Feshbach resonance in combination with 
rf association |9[]. 

In interpreting the experimental results and for future ex- 
tensions, it is essential to develop a detailed understanding of 
the interaction of two particles across the Feshbach resonance, 
taking into account the external confinement of the optical lat- 
tice in a consistent manner. In a seminal paper, Th. Busch et 
al. Hill have analytically solved the problem of two identical 
atoms in a harmonic trap. This model has been compared to 
two-component Fermi gases in an optical lattice at a Feshbach 
resonance Il2lll3ll . 

In this paper, we study the generalized case of two differ- 
ent atoms at an optical lattice site accounting for the anhar- 
monic part of the potential. Both the fact that the two atoms 
feel different trap frequencies and the anharmonicity lead to 
a coupling of center-of-mass and relative motion of the two 
atoms resulting in deviations from the model in Ref. fiUl . We 
model interactions between two cold atoms by a regularized 
6 function type interatomic pseudopotential. We discuss the 
solutions of the uncoupled problem and develop an exact di- 
agonalization approach to the coupling term. In this very gen- 
eral approach, we find significant deviations from the identi- 
cal particle scenario, the strongest effect being observed for 
repulsively interacting atoms with large mass ratios. We fur- 



ther discuss rf association as a method of determining the en- 
ergy spectrum at a heteronuclear Feshbach resonance between 
87 Rb and 40 K in a 3D optical lattice. We compare the the- 
oretical energy spectrum to the experimental results ||9[] and 
discuss methods of precisely determining the Feshbach reso- 
nance center position based on this comparison. Finally we 
calculate the efficiency of rf association and the lifetime of 
heteronuclear 40 K- 87 Rb molecules and find qualitative agree- 
ment with experimental results. 

I. THEORETICAL MODEL 

In order to model interactions within an atom pair, we con- 
sider an interatomic potential given by a regularized S poten- 
tial ITil [l5l [lril [nil . For two atoms of the same kind in an 
isotropic harmonic trap an analytic solution exists II 1 1)1 - Here 
we consider two different atomic species which are confined 
at a single site of a 3D optical lattice. In this case the atoms ex- 
perience different trapping frequencies and the confining po- 
tential has significant anharmonic features. We use the fol- 
lowing Hamiltonian as a starting point 
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Here mi and mi are the masses of the two atoms, 
oji = y/2Vik 2 /m,i are the trapping frequencies obtained us- 
ing the harmonic approximation to the trapping potential, k 
is the wave number and Vi is the depth of the lattice felt by 
atom i, a s is the scattering length, fi = mi ■ m%/M the re- 
duced mass, M = mi + mi the total mass, f = r\ — f< 1 the 
relative position, and r = \r\ — fvj| is the distance between 
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are necessary to accurately approximate the potential of one 
lattice site given by 
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and similar expressions for V^ttice and V^| ice . The first term 
of Eq. (01 gives rise to the harmonic approximation through 
LUi = \j2Vik 2 /nii, and the remainder gives rise to Korr- 

We introduce relative and center-of-mass coordinates 
r = r\ — T2 and R = [m\r\ + m^fi^jM and define the cor- 
responding frequencies: 
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The transformed Hamiltonian consists of three contributions: 
one center-of-mass harmonic oscillator Hamiltonian i? c .m.> 
one term for the relative motion iJ re i containing a harmonic 
oscillator term and the contact interaction, and one last term 
-ffcoupie grouping together all terms which couple relative and 
center-of-mass motion and which arise from the different trap 
frequencies and the anharmonic corrections: 
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Let us first neglect the coupling terms -ff C oupic- In this case, 
the problem separates into relative and center-of-mass motion, 
with the center-of-mass motion given by harmonic oscilla- 
tor wave functions. The Hamiltonian of the relative motion 
H le \ is solved analytically in Ref. Hill and leads to the energy 
structure 
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with a rc i = \Jhj (/i(xJ ro i). Only I = states are considered 
here, since these are the only ones affected by the regularized 
6 potential. The rest of the spectrum consists of I ^ har- 
monic oscillator states with an energy independent of a s . The 
corresponding eigenfunctions are given by 
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FIG. 1: (Color online) Energy eigenvalues of 40 K and 87 Rb as a 
function of scattering length without (black dashed line) and with 
coupling terms (blue solid line) due to anharmonicity and unequal 
trap frequencies in the lattice for parameters: Vnb = 40.5i5 r ,Rb, 
Vk = 0.86Vkb, and A = 1030 nm. The deviation between the ide- 
alized model and the full solution is substantial in particular for the 
upper branch. 



A is a normalization constant which we determine numeri- 
cally, U(a, b; z) are the confluent hypergeometric functions 
and the non-integer indices v are related to the energy by 
E ie] = Huj Iei (2i/ + |). 

The resulting energy spectrum is shown in Fig. [TJ (black 
dashed line) for a center-of-mass energy of 3/2 Huj CAn .. For 
vanishing interaction, the lowest harmonic oscillator state has 
an energy of 3/2(fiw c . m . + hui le \). For large positive values of 
a s , it transforms into repulsively interacting atom pairs with 
a unitary positive "binding energy" of +hoj IC \. In a recent 
experiment with bosonic atoms in an optical lattice, such re- 
pulsively interacting atom pairs served as a starting point to 
create stable repulsively bound pairs |18y . For negative a s , 
the aforementioned state transforms into attractively interact- 
ing atoms. In the unitary limit (a s — » — 00), these atoms ac- 
quire a binding energy of — hx) Te \. When the scattering length 
changes from large and negative to large and positive (as ob- 
served, e. g., at atomic Feshbach resonances), we enter the 
molecule part of the spectrum. In that part of the spectrum, 
the resulting two-body bound state is stable even in the ab- 
sence of the external potential. As a s becomes smaller and 
smaller again from above (a s — ► +0), the size of the molecule 
decreases proportionally to a s , and the binding energy tends 
to —00. 

As soon as we add the coupling term, H coup \ c , this treat- 
ment is no longer valid as center-of-mass and relative mo- 
tion are no longer decoupled. In order to describe this prob- 
lem in a consistent fashion, we have calculated the matrix 
of the complete Hamiltonian (0 using the eigenfunctions of 
H c .m. + H IC \ and numerically obtained energy eigenvalues 
and eigenfunctions for the coupled problem by diagonalizing 
H. 

The anharmonic corrections are treated as follows. Since 
x\ = X + ax and X2 = X — bx with a := m^jM and b := 
mi/M, the cc-dependent part of the anharmonic corrections 
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Corresponding expressions are obtained for the y- and in- 



direction, Vcorr and Vcor'r- In the numerical implementation, 



we have tested for convergence with terms up to eighth order. 
We found that including eighth-order corrections improve the 
accuracy of the calculation by only sw 3 x 10~ 3 fiw re i. 

Our approach leads to the diagonalization of rather small 
Hamiltonian matrices as our main interest is the modi- 
fication of the ground state and the repulsively interact- 
ing pair branch. The whole calculation has been done 
with MATHEMATICA. As basis we have chosen the states 
\N, L, M,n,l,m) with lowest principal quantum numbers 
II := 2N + L + 2n + 1 = 0, 1, n max , where N, L, M 
and n, I, m are the quantum numbers of the eigenfunctions 
of the rotationally symmetric harmonic oscillator of center- 
of-mass and relative motion, respectively. We typically used 
n max = 7 leading to a total number of 258 basis states. We 
have found that adding another level of the uncoupled prob- 
lem to the basis set leads to additional changes in the en- 
ergy smaller than w 10~ 3 fuo IC \. Furthermore, we exploited 
the fact that the total angular momentum L z — H(M + rn) 
of the low-energy eigenfunctions is approximately conserved 
despite the cubic symmetry of the optical lattice. Again, we 
found that including L z j£0 basis states lowers the energy by 
only ps 3 x 10~ 3 fiu/ re i QfJ. 

Fig. [TJ shows the resulting energy spectrum (blue solid 
line) compared to the uncoupled solution (black dashed line), 
calculated for 40 K and 87 Rb with the experimental param- 
eters of Ref. H: V Rh = 40.5£ r , R b, V K = 0.86 V Rh , and 
A = 1030 nm. £7 r ,Rb = fi 2 k 2 /2m R b is the 87 Rb recoil en- 
ergy. In the case of heteronuclear atom pairs it is useful 
to express the lattice depth in units of E r , IC \ — h 2 k 2 /2fx, 
which is the kinetic energy given to a particle with re- 
duced mass /j by a photon of momentum Hk. Then, 
VRb = 40.5 £ r> Rb = 12.6 -E r . r ci. As can be seen from the fig- 
ure, the deviation between the idealized model, where the cou- 
pling term has been neglected, and the full solution is substan- 
tial. The difference is most pronounced in the repulsively in- 
teracting pair branch (0.34 hui IC \ w 20% of the level spacing) 
and becomes smaller as we enter the attractively interacting 
atom part of the spectrum. The molecular branch is relatively 
unaffected by the coupling term H coup i c . This is natural be- 
cause as we approach a — ► +0, the role of the external con- 
finement decreases since the molecule becomes smaller. 

The influence of the coupling term If couple is even stronger 
if we consider molecules with large mass ratios as in the case 
of 6 Li and 133 Cs, see Fig. [2] We have chosen the lattice pa- 
rameters Vu = Vhs = 10 £y,rei an d A = 1 /urn. Here the en- 
ergy of the repulsively interacting atoms is lowered by up to 
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FIG. 2: (Color online) Low-energy spectrum of states with center- 
of-mass energy 3/2 huj c .m. for 6 Li and 133 Cs and lattice parameters 
Vhi = Vcs = 10 h 2 k 2 /2fj, and A = lfim. The energy is much more 



lowered compared to the case of * U K and 
large mass ratio of the 6 Li and 133 Cs atoms. 



Rb. This is due to the 



Table Q] shows the effect of the individual coupling terms, 
Halo '■= /iAcj 2 r • R and V^ rr> on to the energy of several 
atom pairs. The energies have been calculated for repulsively 
interacting atoms at a s = 4 a rc i which is the largest scatter- 
ing length shown in Figs. Q] and [2] All energies of Table U 
are given in units of the level spacing of the relative motion 
hu IC \. Adding the coupling term H^lo contributes up to 62% 
to the total change AE for 6 Li and 133 Cs. The strong influ- 
ence of Haui stems from the large mass ratio which results in 
extremely different trap frequencies uju and wc s - By contrast, 
the energy of 6 Li and 7 Li atoms is nearly not affected by Hau, 
since the trap frequencies are almost equal. 



II. EXPERIMENTAL PROCEDURE 

In the experiment, we have tested our theoretical approach 
by studying the energy spectrum of 40 K and 87 Rb atom 



TABLE I: Influence of the individual coupling terms Ha^ and V CO n, 
onto the total energy of several atom pairs. The energies are given 
in units of hui IC \. All values are calculated at a s — 4a rc i for lat- 
tice depths of Vi = Va = 10 -E r , r ei and a wavelength of A = 1 /im. 
Eo := -Ec.m. + E Ic i is the energy of the uncoupled Hamiltonian. In- 
cluding Hauj into the Hamiltonian reduces the energy by AEaui and 
including Hau> + V CO n reduces the energy further by Ai5 corr . The 
value in brackets is the percentage contribution of the individual cou- 
pling terms to the total change of the energy AE. 
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pairs at a single lattice site of a 3D optical lattice in the 
vicinity of a heteronuclear Feshbach resonance, allowing a s 
to be tuned from strong attractive to repulsive interactions. 
Our experimental procedure for obtaining Fermi-Bose mix- 
tures Il20l J 2III in optical lattices has been discussed previ- 
ously 0, 111]. In our experiment, we obtain a mixture of 
40 K atoms in the \F = 9/2, m F = 9/2) state and 87 Rb in the 
\F = 2, rriF = 2) state by rf-induced sympathetic cooling in a 
magnetic trap. The mixture is transferred into a crossed opti- 
cal dipole trap with final trap frequencies for 87 Rb of 2tt50 Hz. 
In the optical dipole trap, 87 Rb atoms are transferred from 
1 2, 2) to 1 1, 1} by a microwave sweep at 20 G and any remain- 
ing atoms in the upper hyperfine \F = 2,iriF) states are re- 
moved by a resonant light pulse. Next, we transfer 40 K into 
the 1 9/2 , —7/2) state by performing an rf sweep at the same 
magnetic field with almost 100% efficiency. With the mix- 
ture in the 87 Rb|l, 1) ® 40 K|9/2, -7/2) state, we ramp up the 
magnetic field to final field values at the Feshbach resonance 
occurring around 547 G Il22ll23l 12411 . Note that the state which 
we prepare is not Feshbach-resonant at the magnetic field val- 
ues which we study, and that a final transfer of 40 K into the 
|9/2, —9/2) state is necessary to access the resonantly inter- 
acting state. This is precisely the transition which we use to 
measure the energy spectrum as outlined further below. 

We ramp up a 3D optical lattice at a wavelength of 
1030 nm, where the trapping potential for both species is 
related according to VK = 0.86VRb- Due to the differ- 
ent masses of the two species, the trapping frequencies are 
ujk = \/87/40 • 0.86 WRb — 1.4 WRb in the harmonic approx- 
imation. The optical lattice light is derived from a frequency 
stabilized 20W Yb:YAG disc laser with a 50 ms linewidth of 
20 kHz. The lattice is formed by three retroreflected laser 
beams with orthogonal polarizations and a minimum detuning 
of 10 MHz between individual beams. In order to get a maxi- 
mum of lattice sites occupied by one boson and one fermion, 
the best trade-off has been to limit the particle number at this 
stage to a few ten thousand. 

In the optical lattice, we study the energy spectrum (bind- 
ing energy) of two particles at a single lattice site by rf spec- 
troscopy (association) of pairs of one 87 Rb and one 40 K atom 
(see Fig. [3). The idea for the measurement is to drive an rf 
transition between two atomic sublevels one of which is char- 
acterized by the presence of the Feshbach resonance and ex- 
hibits a large variation of the scattering length as a function of 
magnetic field according to l25ll 



a s (B) = a bg ( 1 



AB 

B — Br, 



(11) 



where ab g is the non-resonant background scattering length, 
AB the magnetic field width of the resonance, and Br, the 
resonance center position. The other level involved in the rf 
transition is characterized by a non-resonant scattering length 
independent of magnetic field over the experimentally studied 
field range. We use the 40 K |9/2, -7/2) -> |9/2, -9/2) tran- 
sition where the Feshbach-resonant state is the final |1, 1) ® 
|9/2, -9/2) state. 

A sample spectrum of this transition for the mixture in 
the optical lattice is shown in Fig. [3] The figure shows two 
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FIG. 3: (Color online) rf spectroscopy of 40 K - 87 Rb in a 3D op- 
tical lattice on the 40 K|9/2, -7/2) -> |9/2, -9/2) transition (see 
inset) at a lattice depth of Vkb = 27.5 £ r ,Rb and a magnetic field of 
547.13 G, where the interaction is attractive. The spectrum is plot- 
ted as a function of detuning from the undisturbed atomic resonance 
frequency and clearly shows the large atomic peak at zero detuning. 
The peak at 13.9kHz is due to association of |L 1) <X> 1 9/2, —7/2) 
atom pairs into a bound state (figure from Ref. |9|]). 



peaks; one of them occurs at the frequency corresponding 
to the undisturbed 40 K |9/2, -7/2) -► |9/2, -9/2) Zeeman 
transition frequency at lattice sites occupied by a single 40 K 
fermion. This peak is used for the calibration of the magnetic 
field across the Feshbach resonance using the Breit-Rabi for- 
mula for 40 K and the hyperfine parameters from Ref. l26ll . 
For 57 measurements on 1 1 consecutive days, we find a mean 
deviation from the magnetic field calibration of 2.7 mG at 
magnetic fields around 547 G, corresponding to a field repro- 
ducibility of 5 x 10 ~ 6 . There is an additional uncertainty on 
the absolute value of the magnetic field due to the specified 
reference frequency source accuracy for the rf synthesizer of 
1 x 10~ 5 , resulting in an uncertainty of the measured magnetic 
fields of 12 mG. 

The second peak at a positive detuning of 13.9 kHz is the 
result of interactions between 40 K and 87 Rb at a lattice site 
where one heteronuclear atom pair is present. There are two 
different energy shifts causing the observed separation of the 
peaks: One is the constant, small energy shift of the initial 
|1, 1) ® 1 9/2 , —7/2) state which is independent of B, and 
the important, magnetic field sensitive collisional shift which 
stems from the strong Feshbach-resonant interactions in the 
|1, 1) tg) 1 9/2 , —9/2) final state. In the specific example, the 
negative energy shift (binding energy) of the final state in- 
creases the transition frequency as seen in Fig. [3] In or- 
der to perform spectroscopy on the aforementioned transition, 
we use pulses with a gaussian amplitude envelope (1/e 2 full 
width of 400 /is and total pulse length of 800 /is), resulting 
in an rf 1/e 2 half linewidth of 1.7 kHz. We choose the pulse 
power such as to achieve full transfer on the single atom tran- 
sition, i. e. rf pulse parameters including power are identical 
for all magnetic fields. Not only for rf spectroscopy pulse gen- 
eration, but also for evaporation and state transfer, we have 
used an advanced rf svnthesizer ll52ll allowing precise control 
of frequency, amplitude and phase down to the 5 ns level and 
therefore the synthesis of in principle arbitrary pulse shapes. 
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FIG. 4: (Color online) Experimentally observed energy spectrum to- 
gether with theory without free parameters (black dashed line) and a 
least squares fit for the resonance parameters Bo and AB (red solid 
line). Experimental data from Ref. [9|]. 



HI. EXPERIMENTAL vs THEORETICAL SPECTRUM. 
RESONANCE POSITION 

From rf spectra as in Fig. [3] we can determine the separa- 
tion between the single atom and the two-particle ("molecu- 
lar") peak with high precision (typical uncertainty of 250 Hz) 
and thus extract the binding energy up to a constant offset due 
to non-zero background scattering lengths. At the same time, 
the atomic peak provides us with a precise magnetic field cal- 
ibration as described above. Spectra as in Fig. [3] have been 
recorded for magnetic field values across the whole resonance 
and yield the energy spectrum as a function of magnetic field. 

Fig. [4] shows the measured energy shift across the reso- 
nance at a lattice depth of 40.5_E' r 3b as a function of magnetic 
field. The energy shift is obtained from Fig. [T] by subtract- 
ing the energy of the initial 87 Rb|l J^<g) 40 K|9/2, -7/2) state: 
E S =E - E{a_ 7/2 = -175a ) Q. In addition, Figs. H 
and[TJare connected through Eq. (TTTT) . One branch of the spec- 
trum is characterized by the presence of a "positive" binding 
energy, the repulsively interacting pair branch. In Fig. Q] we 
have seen that this branch continuously transforms into attrac- 
tively interacting atoms as a function of a s . As a function of 
magnetic field, however, and as a result of Eq. (fTTT i. we ob- 
serve this transition as a jump from the left-hand side of Fig.|H 
where the interaction is weak and repulsive, to the right-hand 
side of Fig. |4j where the interaction is weak and attractive. 
Here, we find attractively interacting atoms which decay into 
free atom pairs if the external confinement of the optical lat- 
tice is removed. 

Whereas in Fig.Q] the attractively interacting atoms branch 
and the molecule branch are only asymptotically equal in the 
limit \a s \ — > oo, the singularity on resonance in Eq. ( fTTT i trans- 
forms this into a continuous crossover across the center of the 
resonance position as a function of magnetic field and as seen 
in Fig. [4] These molecules are stable even in the absence of 
the optical lattice potential. 

In order to compare the numerically calculated energy spec- 
trum E(a s ) (blue solid line of Fig.[TJ to the experimental data 
E(B) of Fig. HI we transform the scattering length o s into 



the magnetic field strength B via Eq. ( fTTT i. By using param- 
eters from the literature: a^g = — 185ao, AB = — 3G l28tl 
and Bq = 546.8 G j24ll . we obtain the black dashed curve in 
Fig. [4] As can be seen, the difference between the theoretical 
prediction and the experimental data can be overcome by a 
shift of the black dashed curve along the magnetic field axis. 
We attribute this shift to an insufficient knowledge of the res- 
onance center position Bq. 

We therefore fit our theoretical calculations to the experi- 
mental data in order to improve the estimate for the resonance 
center position Bq. As independent fit parameters we choose 
Bq and AB, while ab g is fixed. The latter parameter cannot 
be determined independently from the measurements due to 
its strong correlations with AB. This is due to the fact that 
in the vicinity of the resonance center position Bq the first 
term of Eq. ( fTTT > is negligible so that only the product a^ g AB 
can be determined precisely from the fit. We therefore set 
abg = — 185ao 128|| and use AB and Bq as free fit parame- 
ters, with the caveat that only the value obtained for Bq is to 
be considered precise. In Fig.|4j the result of the least squares 
fit is displayed as a red solid line. Note that the reliability of 
the fitting procedure sensitively depends on an accurate calcu- 
lation of the energy spectrum E(a s ) which includes an exact 
treatment of the anharmonicity and the different trap frequen- 
cies of the two atoms. 

The least squares fit results in the following values of the 
resonance parameters AB — — 2.92G and Bq = 546.669 G. 
The fit results in an uncertainty of 2 mG on Bq. The value of 
Bq sensitively depends on the scattering length of the initial 
state a_7/ 2 - Assuming an uncertainty of a_7/ 2 of 10% results 
in an uncertainty on Bq of 20 mG. Another possible source of 
systematic uncertainties may be the lattice depth calibration. 
The lattice depth has been calibrated by parametric excitation 
from the first to the third band of the lattice and is estimated to 
have an uncertainty of 5%. Repeating the fit procedure with 
±5% variations on the lattice depth calibration, we obtain a 
corresponding systematic uncertainty on Bq of 7 mG. A third 
source of systematic uncertainties finally stems from the finite 
basis and an imprecise approximation of the lattice site po- 
tential. Here, we included corrections up to eighth order and 
generated basis states of the lowest eight energy levels of the 
uncoupled Hamiltonian. This improved the value of Bq by 
2 mG compared to a calculation with up to sixth order correc- 
tions and basis states of lowest seven energy levels. Adding 
the systematic uncertainty of the magnetic field calibration of 
12mG (see above), we finally obtain 

B = 546.669(24) syst (2) stat G (12) 

under the assumption that the pseudopotential treatment is 
valid in the present experimental situation Hill . 

IV. EFFICIENCY OF rf ASSOCIATION 

In a next step, we have analyzed the transfer efficiency of 
the rf association. The rf association process can be described 
theoretically by a Rabi model: The spin of the 40 K atoms is 
flipped from |9/2, -7/2} to |9/2, -9/2) by applying a radio 
frequency. The atoms are initially in state |1) := 0) and 
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afterwards in state 2} := (0, $/). In the rotating frame and 
by integrating out the spatial degrees of freedom we obtain 
the Hamiltonian 



H I -Auj 
2 \ (^ f ) Ul (t) 



<*i|$/)wi(f) 
+Auj 



(13) 



where Auj := uj — ujq — u>b is the detuning, uj is the radio fre- 
quency, ujq oc B n is proportional to the applied magnetic field, 
u>b is proportional to the binding energy £5 = fajJb, is 
proportional to the amplitude of the oscillating magnetic field 
Bi(t), and is the overlap integral between the ini- 

tial and final motional wave functions. As can be seen, the 
off-diagonal elements of Hamiltonian (l3[ are not only pro- 
portional to the rf amplitude u>i(t), but also to the overlap 
integral ($;| <£>/}. Therefore, the transfer probability between 
states |1) and 2} corresponds to Rabi flopping with a Rabi 
frequency reduced by the overlap integral of $i and $ t com- 
pared to the pure atomic transition. Exactly on the molecu- 
lar resonance, we have uj = uj + w& (— > Auj = 0). The on- 
resonant result for the theoretical transfer probability (effi- 
ciency) is thus given by 



Px^ 2 = sin 2 



(14) 



which is unity for a transfer between atomic states, where 
$i = when setting the area under the uj\(t) curve to 

f Q u)i(t')db' = it. For transfer into the molecular state, the 
theoretical probability decreases as a function of the wave 
function overlap integral since the molecular final orbital wave 
function becomes more and more dissimilar from the initial 
two-body atomic wave function. 

In the experiment, the molecules were associated using rf 
pulses designed to induce a it pulse for the non-interacting 
atoms: J Q u>i(t') dt' = tt. This n pulse has been kept fixed 
over the entire range of magnetic field values investigated. 
The experimental association efficiency is determined from 
the height of the molecular peak (see Fig. 3) as a function of 
magnetic field for constant pulse parameters and uj — ujq + ujb 
(— > Auj = 0) as in the theory above. 

Figs. |3 a) and [2b) show a comparison between the con- 
version efficiency as extracted from the experimental data and 
the theoretical estimate from equation (TBI . Theory and ex- 
periment show the general trend of dropping association ef- 
ficiency with increasing binding energy when the initial and 
final wave functions become more and more dissimilar. In 
this context, we define the experimental conversion efficiency 
as the ratio of the number of molecules created and the initial 
lattice sites which are occupied by exactly one K and one Rb 
atom. Note that only on these lattice sites molecules can be 
created. For the comparison of experimental and theoretical 
transfer efficiency, the experimental data have been scaled by 
a global factor to reproduce a conversion efficiency of 1 far 
off the Feshbach resonance where initial and final two-body 
wave function are equal. This is necessary, because the initial 
lattice sites occupied by one K and one Rb atom have not been 
determined experimentally. 
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FIG. 5: (Color online) Transfer efficiency of rf association as ob- 
served in the experiment and estimated from a Rabi model, both for 
(a) attractively interacting atoms and molecules, and for (b) repul- 
sively interacting pairs. The experimental data contain a global factor 
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While the experiments presented here were performed at 
constant rf pulse parameters, it should be possible from the 
above arguments to increase either pulse power or duration or 
both of the rf pulse to account for the reduced wave function 
overlap and thereby always obtain an efficiency of 1 . In par- 
ticular, it should be possible to drive Rabi oscillations between 
atoms and molecules in a very similar way as recently demon- 
strated 1 29]. The comparison also indicates that in the case 
of association efficiency a full quantitative agreement might 
require some more sophisticated treatment of the association 
process. This is in contrast to the analysis of binding energies 
and lifetimes (see below), where the good quantitiative agree- 
ment shows that here the 6 interaction approximation captures 
the essential physics. Testing the Rabi oscillation hypothesis 
for molecules with rf might provide further insight. 

V. LIFETIME 

Molecule formation at atomic Feshbach resonances re- 
sults in dimers which are very weakly bound and may ex- 
hibit strong inelastic collisional losses. Experiments with 
molecules created from bosonic atoms showed very small 
lifetimes. As a result, these molecules can be brought into 
the quantum degenerate regime 13011 . but thermal equilibrium 
is generally difficult to achieve for molecules created from 
bosonic atoms because of the short lifetime. 
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FIG. 6: (Color online) Sketch of expected lattice occupation. The 
arrows illustrate tunneling of remaining fermionic 40 K atoms to the 
"molecular" shell where they can undergo inelastic three-body colli- 
sions with a 4() K- 87 Rb molecule. 



In experiments with molecule creation from two- 
component Fermi gases OH], inelastic molecule-molecule 
and molecule-atom collisions are suppressed by the Pauli 
exclusion principle IL321 l33ll . resulting in remarkably long 
lifetimes between approximately 100 ms and even sec- 
onds, allowing Bose-Einstein condensation of Feshbach 
molecule s B3411 a nd the observation of BCS-BEC crossover 
physics i35i [M |H HI & The lifetime limitation for 
molecules created from bosonic atoms has been overcome by 
creating molecules in 3D optical lattices, where molecules 
are created at a single lattice site and isolated from inelastic 
collisions with residual atoms or other molecules l4(ill . 

For molecules created from Fermi-Bose mixtures, the sit- 
uation is a little bit more complicated. As far as collisions 
between molecules are concerned, the fermionic character 
of the molecule should become more evident the deeper the 
molecule is bound, thus resulting in suppression of colli- 
sions [41]. 

As far as collisions with residual atoms are concerned, we 
expect inelastic collisions with fermionic atoms in the same 
spin state as the fermionic component of the molecule, i. e. in 
the 1 9/2 , —9/2) state, to be suppressed due to the Pauli exclu- 
sion principle close to the resonance, when the "atomic" char- 
acter of the molecule's constituents is still significant ll32l[33ll . 
For collisions with bosonic atoms and fermionic atoms in a 
different internal state, we do not expect any Pauli-blocking 
enhanced lifetime, since the residual atom can in principle 
come arbitrarily near to the molecule's constituents. 

In our situation, where the molecules are created through 
rf association, residual fermionic atoms remain in a different 
spin state, either in |9/2, -7/2) or |9/2, -5/2) (for the latter 
case, and for a description of the experimental procedure, see 
Ref. H). These residual fermionic atoms as well as the re- 
maining bosonic atoms may therefore limit the stability of the 
molecular sample. 

Molecule creation in the optical lattice introduces a second 
aspect concerning the lifetime: lattice occupation and tunnel- 
ing probabilities. In Fig. [6] we have sketched the expected 
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FIG. 7: (Color online) Lifetime of heteronuclear 40 K- 87 Rb 
molecules in the optical lattice. The Lifetime is limited due to resid- 
ual atoms which can tunnel to lattice sites with molecules and pro- 
voke inelastic three-body loss. The theoretical prediction uses the 
pseudopotential wave function and contains a global factor which 
was adjusted to the experimental data of Ref. OJ. 

occupation in the optical lattice. Prior to molecule creation, 
we expect slightly less than unity filling for the fermionic 
component. As far as the bosons are concerned, we expect 
a central occupation number between 3 and 5, surrounded by 
shells of decreasing occupation number. In the rf associa- 
tion process, molecules are only created in the shell where 
we have one fermion and one boson per lattice site. In the 
outermost region of the lattice, we have lattice sites with only 
one fermion which are responsible for the "atomic" peak in 
the rf spectroscopy signal. After the rf association process in 
the "molecular" shell, bosons from neighboring sites as well 
as fermions remaining in a different spin state can tunnel to 
the "molecular" shell and provoke inelastic three-body loss. 
In our experimental situation, this is more probable for the 
remaining fermionic atoms, since they are lighter and have a 
smaller tunneling time (10 to 20 ms for the lattice depths dis- 
cussed here). For the highest binding energies observed in 
the experiment, we find a limiting lifetime of 10 to 20 ms as 
seen in Fig. [7] which is consistent with the assumption that 
in this case, three-body loss is highly probable once tunneling 
of a distinguishable residual fermion has occured. Still, for 
the more weakly bound molecules and in particular for attrac- 
tively interacting atoms, we observe high lifetimes of 120 ms, 
raising the question of the magnetic field dependence of the 
lifetime. 

We can understand this magnetic field dependence using 
the pseudopotential model by introducing a product wave 
function for the combined wave function of the resonantly in- 
teracting atom pair and a residual fermionic atom after tunnel- 
ing to a molecular site. We write this three-body wave func- 
tion as 

$(r, R, r 3 ) := $ mol (r, #)$ 3 (f 3 ) (15) 

where $ mo i is the result of the pseudopotential calculation for 
the molecule and $3 is the ground-state wave function of the 
residual atom at the same lattice site. Note that this treatment 
assumes weak interactions between the residual atom and the 
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molecule (the interaction between the residual atom and the 
molecule's constituents is on the order of the background 
scattering length). From solution (l5[ of the pseudopotential 
model, the dependence of the loss rate on the scattering length 
can be obtained close to the resonance 113211 : the loss rate Y is 
proportional to the probability V of finding the three atoms 
within a small sphere of radius a, where they can undergo 
three-body recombination. This probability is expected to be- 
come larger for more deeply bound molecules, since two of 
the three atoms are already at a close distance. Up to a global 
factor, V is independent of the value chosen for er, provided 
(7 a ro i, and also a <C a s in the molecule regime. More 
quantitatively, we calculate this probability according to 

V= I drdRdr 3 \3>{r,R,r 3 )\ 2 . (16) 

J \r\«r 1 1 

l?3-«l<" 

The magnetic field dependence of the loss rate is thus 
given through Y oc P, and the lifetime is proportional to 
1/r. By using the wave functions (O and the relation 

\ A\ 2 oc a 2 dE le i/da s Hill , we obtain 



2ftcd tel ^ i> V\ 2tm tB i T i> 

where C is independent of a s , and ij){x) = Y'{x)/Y(x) is the 
digamma function. This result agrees with a numerical inte- 
gration of Eq. ( fTol l using the eigenfunctions of the complete 
Hamiltonian (0. 

The lifetime obtained from the calculation is shown in 
Fig. [7] as a red solid line, scaled by a global factor to allow 
comparison to the experiment. As can be seen, the theoretical 
prediction explains the magnetic field dependence of the life- 
time rather well. From an experimental point of view, we can 
therefore expect that removal of the remaining atoms using a 
resonant light pulse will significantly increase the lifetime of 
the molecules in the optical lattice. 



CONCLUSIONS 

To summarize, we have developed a pseudopotential ap- 
proach to the scattering of unequally trapped atoms at a single 
site of an optical lattice including terms which couple center- 
of-mass and relative motion. We have compared the energy 
spectrum to experimental results for 40 K and 87 Rb atoms in- 
teracting at a heteronuclear Feshbach resonance in a 3D op- 
tical lattice. Within the pseudopotential model we have pre- 
cisely determined the center position of the Feshbach reso- 
nance based on this comparison. The pseudopotential ap- 
proach also allows us to understand the efficiency of rf as- 
sociation used to experimentally determine the energy spec- 
trum, as well as the dependence of the molecular lifetime on 
magnetic field. The model developed in this paper enables a 
broad understanding of heteronuclear atom pairs in an opti- 
cal lattice. We are aware of possible limitations of the pseu- 
dopotential model. It might be an interesting option to ex- 
tend the method described here to energy dependent pseu- 
dopotential s lil 51 BE S3l or multi-channel mod- 
els 11431 l48l 14911 . Finally, we note that the present rf associ- 
ation technique could be used to study the three-body prob- 
lem at a triply occupied lattice site lf50l l5lll . An advantage of 
this method with respect to the adiabatic magnetic field sweep 
proposed in 1 50] is that it is less sensitive to three-body losses, 
which are particularly important for Efimovian states [51]. 

ACKNOWLEDGMENTS 

The authors would like to thank K. Rzazewski for his con- 
tributions in the initial discussion of possible extensions to 
Ref. fiUl . We acknowledge discussions with Y. Castin, R Juli- 
enne, Th. Kohler, R Naidon, D. Petrov, and L. Pricoupenko on 
the pseudopotential models and on loss rates. We would like 
to thank A. Simoni for providing us with the closed channel 
magnetic moment at the Feshbach resonance and the scatter- 
ing length in the initial state of the rf association. We ac- 
knowledge financial support by the Deutsche Forschungsge- 
meinschaft (SPP 1116). 



[1] L. Mathey, D.-W. Wang, W. Hofstetter, M. D. Lukin, and Eu- 
gene Dernier, Phys. Rev. Lett. 93, 120404 (2004); E. Pazy and 
A. Vardi, Phys. Rev. A 72, 033609 (2005); D.-W. Wang, M. D. 
Lukin, and E. Demler, Phys. Rev. A 72, 051604(R) (2005). 

[2] L. Viverit, Phys. Rev. A 66, 023605 (2002). 

[3] M. J. Bijlsma, B. A. Heringa, and H. T. C. Stoof, Phys. Rev. A 
61, 053601 (2000). 

[4] D. V. Efremov and L. Viverit, Phys. Rev. B 65, 134519 (2002). 

[5] F. Illuminati and A. Albus, Phys. Rev. Lett. 93, 090406 (2004). 

[6] H. P. Biichler and G. Blatter, Phys. Rev. Lett. 91, 130404 
(2003). 

[7] S. Ospelkaus, C. Ospelkaus, O. Wille, M. Succo, P. Ernst, K. 
Sengstock, and K. Bongs, Phys. Rev. Lett. 96, 180403 (2006). 

[8] Kenneth Giinter, Thilo Stoferle, Henning Moritz, Michael 
K6hl, and Tilman Esslinger, Phys. Rev. Lett. 96, 180402 (2006). 

[9] C. Ospelkaus, S. Ospelkaus, L. Humbert, P. Ernst, K. Seng- 
stock, and K. Bongs, Phys. Rev. Lett. 97, 120402 (2006). 
[10] S. B. Papp and C. E. Wieman, Phys. Rev. Lett. 97, 180404 
(2006). 



[11] Th. Busch, B. G. Englert, K. Rzazewski, and M. Wilkens, 
Found. Phys. 28, 549 (1998). 

[12] T. Stoferle, H. Moritz, K. Giinter, M. Kohl, and T. Esslinger, 
Phys. Rev. Lett. 96, 030401 (2006). 

[13] M. Kohl, K. Giinter, T. Stoferle, H. Moritz, and T. Esslinger, J. 
Phys. B, S47-S56, 39 (2006). 

[14] We expect the pseudopotential model to be fairly accurate for 
the experimental parameters of [ 9]. Indeed, even for large scat- 
tering lengths, this model is expected to become exact in the 
zero-range limit |15], that is when l/fet yp 3> max(/36, |r2|)- 
Here, k typ is the typical wavenumber of the relative motion of 
the two atoms, /?e = (2/xCe/ft 2 ) 1//4 is the van der Waals length 
scale, and r° := -h 2 /(na hg ABdE reB /dB), 8E ICB /dB be- 
ing the magnetic moment of the closed channel with respect 
to the two-atom open channel. For K-Rb in their ground state, 
fa = 7.6 nm (H. Using 8E ICB /dB = fc s 144^K/G (H, we 
get r° = -4.6 nm. It remains to estimate ktyp. In the molecule 
regime, we have l/fc typ ~ a s > 47 nm. In the other regimes 
(attractively and repulsively interacting atoms), we have 



9 



1/fctyp ~ a rci- In the harmonic approximation for the experi- 
mental lattice depth, a ro i = 103 nm. Thus, the above inequality 
is fairly well satisfied. 

[15] Y. Castin, cond-mat/0612613 (to be published in the Proceed- 
ings of the Enrico Fermi School on Fermi gases). 

[16] A. Derevianko, J. F. Babb, and A. Dalgarno, Phys. Rev. A 63, 
052704 (2001). 

[17] A. Simoni, private communication. 

[18] K. Winkler, G. Thalhammer, F. Lang, R. Grimm, J. Hecker 
Denschlag, A. J. Daley, A. Kantian, H. P. Biichler, P. Zoller, 
Nature 441, 853 (2006); A. J. Daley, A. Kantian, H. P. Biichler, 
P. Zoller, K. Winkler, G. Thalhammer, F. Lang, R. Grimm, J. 
Hecker Denschlag, Proceedings of ICAP-2006 (Innsbruck) and 
cond-mat/0608721v2. 

[19] We note that we neglect many low-energy basis states 
|JV,L,0,0,0,0> with 2iV + £>n max when the scatter- 
ing length is small and positive a s 0. These molecule 
states with highly excited center-of-mass energy have nev- 
ertheless low total energy since the energy of the rela- 
tive motion is large and negative. We have included the 
states |JV,L,0,0,0,0) with 2N + L < II max to allow for 
the flattening of the center-of-mass wave function due to 
the anharmonicity of the lattice site potential. But we have 
neglected the states with 2JV + L > II max since they are 
only weakly coupled to the states |0, 0, 0, 0, 0, 0) (molecule) 
and |0, 0, 0, n — 1, 0, 0) (repulsively interacting atoms of 
the decoupled problem). These matrix elements become 
even zero if 2JV + L > p (p — 1,...,8) since in this case 
the matrix elements (0, 0, 0, 0, 0, 0\X p \N, L, 0, 0, 0, 0) and 
(0, 0, 0, 1, 0, 0\X p \N, L, 0, 0, 0, 0) are exactly zero. 

[20] C. Ospelkaus, S. Ospelkaus, K. Sengstock, and K. Bongs, Phys. 
Rev. Lett. 96, 020401 (2006). 

[21] S. Ospelkaus, C. Ospelkaus, R. Dinter, J. Fuchs, M. Nakat, K. 
Sengstock, K. Bongs, J. Mod. Opt. 54, 661 (2007). 

[22] S. Inouye, J. Goldwin, M. L. Olsen, C. Ticknor, J. L. Bohn, and 
D. S. Jin, Phys. Rev. Lett. 93, 183201 (2004). 

[23] Francesca Ferlaino, Chiara D'Errico, Giacomo Roati, Matteo 
Zaccanti, Massimo Inguscio, and Giovanni Modugno, Andrea 
Simoni, Phys. Rev. A 73, 040702(R) (2006). 

[24] S. Ospelkaus, C. Ospelkaus, L. Humbert, K. Sengstock, and K. 
Bongs, Phys. Rev. Lett. 97, 120403 (2006). 

[25] T. Kohler, K. Goral, and P. S. Julienne, Rev. Mod. Phys. 78, 
1311 (2006). 

[26] E. Arimondo, M. Inguscio, and P. Violino, Rev. Mod. Phys. 49, 
31 (1977). 

[27] For the scattering length in the initial 
S7 Rb|l, 1) ® 40 K|9/2, -7/2) state in the considered mag- 
netic field range 544G < B < 549G we take the value 
— 175ao (A. Simoni, private communication). 

[28] M. Zaccanti, C. D'Errico, F. Ferlaino, G. Roati, M. Inguscio, 
and G. Modugno, Phys. Rev. A 74, 041605(R) (2006). 

[29] N. Syassen, D. M. Bauer, M. Lettner, D. Dietze, T. Volz, S. Diirr 
and G. Rempe, Phys. Rev. Lett. 99, 033201 (2007). 

[30] K. Xu, T. Mukaiyama, J. R. Abo-Shaeer, J. K. Chin, D. E. 
Miller, and W. Ketterle, Phys. Rev. Lett. 91, 210402 (2003); 
Jens Herbig, Tobias Kraemer, Michael Mark, Tino Weber, 
Cheng Chin, Hanns-Christoph Nagerl, and Rudolf Grimm, Sci- 
ence 301, 1510 (2003). 



[31] C. A. Regal, C. Ticknor, J. L. Bohn and D. S. Jin, Nature 424, 
47 (2003); J. Cubizolles, T. Bourdel, S. J. J. M. F. Kokkel- 
mans, G. V. Shlyapnikov, and C. Salomon, Phys. Rev. Lett. 91, 
240401 (2003); K. E. Strecker, G. B. Partridge and R. G. Hulet, 
Phys. Rev. Lett. 91, 080406 (2003); S. Jochim, M. Bartenstein, 
A. Altmeyer, G. Hendl, C. Chin, J. Hecker Denschlag, and R. 
Grimm, Phys. Rev. Lett. 91, 240402 (2003). 

[32] D. S. Petrov, C. Salomon, and G. V. Shlyapnikov, Phys. Rev. 
Lett. 93, 090404 (2004); Phys. Rev. A 71, 012708 (2005). 

[33] D. S. Petrov, C. Salomon, and G. V. Shlyapnikov, J. Phys. B 38, 
S645-S660 (2005). 

[34] M. Greiner, C. A. Regal and D. S. Jin, Nature 426, 537 (2003); 
M. W. Zwierlein, C. A. Stan, C. H. Schunck, S. M. F. Rau- 
pach, S. Gupta, Z. Hadzibabic, and W. Ketterle, Phys. Rev. Lett. 
91, 250401 (2003); S. Jochim, M. Bartenstein, A. Altmeyer, G. 
Hendl, S. Riedl, C. Chin, J. Hecker Denschlag, and R. Grimm, 
Science 302, 2101 (2003). 

[35] C. A. Regal, M. Greiner, and D. S. Jin, Phys. Rev. Lett. 92, 
040403 (2004). 

[36] M. W. Zwierlein, C. A. Stan, C. H. Schunck, S. M. F. Raupach, 
A. J. Kerman, and W. Ketterle, Phys. Rev. Lett. 92, 120403 
(2004). 

[37] M. Bartenstein, A. Altmeyer, S. Riedl, S. Jochim, C. Chin, J. 
Hecker Denschlag, and R. Grimm, Phys. Rev. Lett. 92, 120401 
(2004). 

[38] T. Bourdel, L. Khaykovich, J. Cubizolles, J. Zhang, F. Chevy, 

M. Teichmann, L. Tarruell, S. J. J. M. F. Kokkelmans, and C. 

Salomon, Phys. Rev. Lett. 93, 050401 (2004). 
[39] J. Kinast, S. L. Hemmer, M. E. Gehm, A. Turlapov, and J. E. 

Thomas, Phys. Rev. Lett. 92, 150402 (2004). 
[40] G. Thalhammer, K. Winkler, F. Lang, S. Schmid, R. Grimm, 

and J. Hecker Denschlag, Phys. Rev. Lett. 96, 050402 (2006); 

T. Volz, N. Syassen, D. M. Bauer, E. Hansis, S. Diirr, G. Rempe, 

Nature Physics 2, 692 (2006). 
[41] C. A. Stan, M. W. Zwierlein, C. H. Schunck, S. M. F. Raupach, 

and W. Ketterle, Phys. Rev. Lett. 93, 143001 (2004). 
[42] D. Blume and C. H. Greene, Phys. Rev. A 65, 043613 (2002). 
[43] E. L. Bolda, E. Tiesinga, and P. S. Julienne, Phys. Rev. A 66, 

013403 (2002). 

[44] Z. Idziaszek and T. Calarco, Phys. Rev. A 74, 022712 (2006). 

[45] D. S. Petrov, Phys. Rev. Lett. 93, 143201 (2004). 

[46] P. Naidon, E. Tiesinga, W. Mitchell, and P. Julienne, New J. 

Phys. 9, 19 (2007); and references therein. 
[47] For an equivalent energy-independent formulation of such a 

model, see: L. Pricoupenko, Phys. Rev. A 73, 012701 (2006). 
[48] E. Tiesinga, C. J. Williams, F. H. Mies, and P. S. Julienne, Phys. 

Rev. A 61, 063416 (2000). 
[49] D. B. M. Dickerscheid and H. T. C. Stoof, Phys. Rev. A 72 

053625 (2005); K. B. Gubbels, D. B. M. Dickerscheid, and H. 

T. C. Stoof, New J. Phys. 8, 151 (2006); and references therein. 
[50] Martin Stoll, Thorsten Kohler, Phys. Rev. A 72, 022714 (2005). 
[51] F. Werner and Y. Castin, Phys. Rev. Lett. 97, 150401 (2006). 
[52] VFG-150, a FPGA-driven fast Versatile-Frequency-Generator, 

controlled via an USB2.0 interface and developed by Th. Han- 

nemann in the group of C. Wunderlich, University of Siegen, 

Germany. 



